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ABSTRACT 

Using deep infrared observations conducted with the MOIRCS imager on the Subaru Telescope in the northern 
GOODS field combined with pubhc surveys in GOODS-S, we investigate the dependence on stellar mass, M*, and 
galaxy type of the close pair fraction {5h~^ kpc < rgep < 20h~^ kpc) and implied merger rate. In terms of combined 
depth and survey area, our publicly available mass-limited sample represents a significant improvement over earlier 
infrared surveys used for this purpose. In common with some recent studies we find that the fraction of paired 
systems that could result in major mergers is low (~4%) and does not increase significantly with redshift to z ~ 1.2, 
with cx (1 + z)^-^*^-^. Our key finding is that massive galaxies with Af, > lO^^M© are more likely to host merging 
companions than less massive systems (M* ~ IO^'^Mq). We find evidence for a higher pair fraction for red, spheroidal 
hosts compared to blue, late-type systems, in line with expectations based on clustering at small scales. So-called "dry" 
mergers between early-type galaxies devoid of star formation represent nearly 50% of close pairs with M, > 3 x lO^^M© 
at z ~ 0.5, but less than 30% at z ~ 1. This result can be explained by the increasing abundance of red, early-type 
galaxies at these masses. We compare the volumetric merger rat e of galaxies with dif ferent masses to mass-dependent 
trends in galaxy evolution. Our results reaffirm the conclusion of iBundv et al.l (|2007f ) that major mergers do not fully 
account for the formation of spheroidal galaxies since z ~ 1. In terms of mass assembly, major mergers contribute 
little to galaxy growth below ~ 3 x 1O^°M0 but play a more significant role among galaxies with M, > lO^^M© 
~30% of which have undergone mostly dry mergers over the observed redshift range. Overall, the relatively rapid and 
recent coalescence of high mass galaxies mirrors the expected hierarchical growth of halos and is consistent with recent 
model predictions, even if the top-down suppression of star formation and morphological evolution (i.e., "downsizing") 
involves additional physical processes. 

Subject headings: galaxies: evolution — galaxies: formation — galaxies: interactions 
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1. INTRODUCTION 

Recent surveys have challenged the traditional view 
that merging between dark matter halos governs the 
mass assembly history and evolution of the galax- 
ies hosted by these halos. The surprising abun- 
dance o^_estaWislied massi v e galaxies at z 2 
(e.g., iGlazebrook et all |2004 ICimatti et all |2004[ ) has 
been followed by increasing evidence for an early 
completion of star formation in the most massive 
galaxies, followed by continued activity in lower-mass 
systems — a ph enomeno n now termed "downsizing" (e.g 



Juneau et al][2 005: Tre u et al.ll2005t iBundv et al.| 



2007 



Borch et a l. 2006; Cim atti et all l2006t iBell et al.l 
Cowie fc Ba rger 200_8). These observations reveal top 
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down evolutionary patterns that stand in contrast to the 
expected hierarchical, or bottom-up, nature of CDM halo 
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assembly. 

While semi-analytic models employing various energy 
feedback prescriptions have moved closer to reproduc- 
ing top-down behavior within the hierarchical framework 
(ICroton et al. 2006; Bower ct al." '2006'; 'Dc L ucia et al.l 
120061 ICattane o et al. 2008: Stringer ct al. 200^, it has 
not been observationally possible to verify the CDM pre- 
diction that via merging, more massive galaxies are as- 
sembled at later times. In principle, the signal is en- 
coded in the evolving stellar mass function (MF) which 
would be expected to show a rising number density at 
the high-mass end as a function of cosmic time. Current 
surveys, however, are too much affected by cosmic vari- 
ance (or, more accurately, sample variance) un certainties 
to confirm this trend (see IStringer et al]|2008f) . It is also 
difficult to separate evolution in the MF resulting from 
galax y mergers from that ass ociated with star formation 
(SF) (|Drorv fc Alvare j|2008D . The clear alternative is to 
study galaxy mergers directly and to quantify their role 
in driving mass assembly by determining how the galaxy 
merger rate depends on mass. 

In addition to driving mass assembly, major mergers 
(defined here to be those involving components with a 
mass ratio greater than 1/4) have also been proposed as 
the key mechanism that transforms disk-like galaxies into 
spheroidals (Toomre 1977). It may also s hut off star for- 
matio n via triggered AGN feedback (see I Hopkins et al.l 
|2007() . Mergers may not only move galaxies onto the 
"red sequence," but so-called "dry mergers" between red- 
sequence systems have been invoked to explain the in- 
creasing contribution at > lO^^M© to t he global stel- 
lar mass density from red galaxies, (e.g.. Ivan Dokkuml 



2 



Bundy, Fukugita, Ellis, Targett, Belli, Kodama 



[200llFaber eraII[2007l : lBell et al.l[200l . 

The fundamental question is whether the major merger 
rate is high enough to explain the mass-dependent in- 
crease in the numbers of spheroidal and red-sequence 

galaxies or whether other mechanisms for bu ilding 

spheroidal systems are needed. .Bundv et al.l ()2007i) ap- 
pealed to numerical simulations to argue that the pre- 
dicted major merger rate among dark matter halos is 
too low. A sim ilar conclusion was also reached by 
iGenel et al.l((2008h . It is now imperative to compare the 
formation rate of spheroidals to observations capable of 
identifying major mergers and resolving their frequency 
as a function of mass. 

Previous observational attempts to measure the galaxy 
merger rate since z ~ 1 have faced a number of chal- 
lenges and sometimes disagree in their conclusions. The 
derived rate depends on the assumed efficiency and 
merger timescale, but more importantly, merger defi- 
nitions vary among authors. Additionally, it has been 
difficult to distinguish major from minor mergers. At 
z < 1, the merger rate has been estimated based 
on th e occurrence of morphologically d i sturbe d systems 
(e.g., IConselice et"all l2003t iLotz et al.l f2008bl ) and the 
frequency of optical pairs, either defined with respect 
to their re lative velocity difference fe.g.. [ Patton et al.l 
I2002t iLin et al, 2Q04, 2008; do Ra vel et al II2008D or cor- 
rected for chance projections (e.g., Le Fevre et al.l[2000l : 
iBell et al.|[2006al: iKartaltepe et al.ll2007D . Parameteriz- 
ing the redshift dependence as cx (1 -|- z)"*, the range of 
reported values is virtually unconstrained, with m — 0- 
4. Part of the problem may be the use of optical samples 
which we argued in previo us work based on a near-IR 
sample (|Bundve^l.ll2004D. may be biased by triggered 
SF (also see iBerrier et al.ll2006( ). 

Although great care is needed in making comparisons 
with these methods, each of which defines a "merger" in 
a different way, it is now clear that there is considerable 
uncertainty in the literature regarding both the rate of 
merging and whether it rises substantially with redshift. 
Our goal in this paper is to overcome some of the prob- 
lems faced in previous work by using a mass-limited sam- 
ple of galaxy pairs drawn from the Gre at Observatories 
Origi ns Deep Survey fields (GOODS, iGiavalisco et al.l 
120041 ). Our Ks'Selected sample is more than 15 times 
larger than that discussed in iBundv et al.l (|2004[ ). The 
enlarged sample allows us to identify pairs using the 
physically motivated definition of a major merger — as 
determined by the inferred mass ratio — and therefore 
enables us to measure the impact of mergers on mass 
assembly and make direct comparisons to evolutionary 
trends such as the formation rate of spheroidal galaxies. 
With our large sample size, we can study the merger rate 
as a function of mass, a key tool for testing the late as- 
sembly times predicted for massive galaxies. We can also 
distinguish dry and "wet" mergers (those involving sig- 
nificant cold gas reservoirs) as well as the types of galax- 
ies most likely to host close companions. This project 
combines publicly available data in GOODS-S with new 
observations obtained in GOODS-N with the recently- 
commissioned Multi -Object Infrared Cam era and Spec- 
trograph fMOIRCS. Ilchikawa et aLl |2006f) panoramic in- 
frared imager on the Subaru 8.2m telescope. 

Using this data set, we address two important ques- 
tions in this paper. The first is whether the rate of ma- 



jor mergers since z ~ 1 is sufficient to explain the for- 
mation of new spheroidal and red-sequence galaxies over 
this redshift range. The second goal is to quantify the 
role of mergers in galaxy growth and determine whether 
the mass dependence of the merger rate is consistent with 
expectations based on hierarchical assembly. 

A plan of the paper follows. In §[5]we introduce the new 
MOIRCS dataset and discuss its image processing and 
simulations undertaken to determine its limiting magni- 
tude. We also present our analysis of complementary 
near-IR data that is publicly available and was obtained 
in GOODS-S using the Infrared Spectrometer and Array 
Camera (ISAAC) on the Very Large Telescope (VLT). 
The resulting i^s-band catalogs are correlated with a 
number of public redshift and imaging surveys, and the 
matched catalogs will be made available online. In §[3l we 
discuss our estimates of stellar masses, rest-frame colors, 
and morphologies. We present our methods for measur- 
ing pair fractions and comoving number densities in § 
m The results are described in § [5l while in § |6] we dis- 
cuss the inferred merger rates and implications for both 
mass assembly and morphological/color transformations. 
A summary is provided in § [T] 

Throughout this paper, we use the AB magnitude sys- 
tem and adopt a standard cosmology with Hq=70 hjQ 
km s-i Mpc^i, nM=0.3 and f^A^O.T. 

2. DATA AND CATALOGS 

Our study of galaxy pairs makes use of the GOODS 
fields, which not only provide deep HST imaging in four 
bands, but have also been targeted by many foUowup 
surveys at a variety of wavelengths. We begin this section 
with a description of the acquisition and reduction of our 
new MOIRCS data that delivers much-needed deep Kg- 
band coverage across the entire GOODS-N. We then turn 
to publicly available data, presenting first our analysis 
of near-IR imaging in GOODS-S obtained with ISAAC. 
These data are very similar in depth and resolution to 
the MOIRCS data, and together the two datasets provide 
the -ftTs-detected source lists that form the basis of the 
samples used in this paper. 

We then discuss how we cross reference our ifg-selected 
catalogs to the many imaging and spectroscopic datasets 
that are available in both GOODS fields. Because a key 
part of our merger rate analysis uses redshift informa- 
tion to confirm potential galaxy pairs, it is important to 
make use of as many spectroscopic redshifts as possible. 
For this reason, we match our sample to the most up-to- 
date versions of available surveys including many of the 
recently completed ESO surveys in GOODS-S^. Finally, 
where previously measured spectroscopic or photomet- 
ric redshifts are not available, we describe the additional 
photometric redshifts that are required to supplement 
our sample. Recognizing the value of ifg-selected multi- 
wavelength compilations in GOODS, the final matched 
catalogs will become publicly available at this website in 
early 2009: astro.berkeley.edu/--kbundy/KGODDS/ 

2.1. MOIRCS Near-IR Imaging 

We obtained Kg-hand observations using MOIRCS on 
the Subaru Telescope during two nights in April 2006. 

^ In GOODS-S, our catalog serves as a partial update to the 
MUSIC compilation of IGrazian et al.l ||2003 '). 
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Each MOIRCS pointing is imaged onto two overlapping 
detectors, giving a total field of view that is roughly 4' 
by 7'. Our observations consisted of 8 tiled pointings 
arranged with some overlap in order to cover the en- 
tire GOODS-N region. At each pointing, we executed 
three to four 9-point dither patterns (dithering by 15"), 
coadding 4 exposures at each position. Individual expo- 
sure times varied from 50 to 55 seconds. The total inte- 
gration time was typically one hour across the field with 
average seeing of 0'.'5. The MOIRCS A'^-band 80% com- 
pleteness depth varies across the GOODS-N field from 
Ks = 23.5 to Kg = 24.0 with an average of Ks ~ 23.8. 

We reduced the data using the MCSRED IRAF pack- 
age written by Ichi Tanaka. Because MOIRCS is a rela- 
tively new instrument, we report in the Appendix on our 
modifications to MCSRED and solutions to other prob- 
lems we encountered during the data reduction. Once 
the final images were reduced and combined we regis- 
tered them to the GOODS-N ACS astrometry by com- 
paring bright stars detected in both the MOIRCS and 
ACS data. Our Subaru observations were not taken in 
completely photometric conditions, and so photometric 
calibration was carried out by comparing to stars in the 
shallower ( Kg 5, 22.7) Paloma r Kg-hanA observations 
presented inl Bundv et a l.' ("2005'). 

We used SExtractor {Bertin fc Arnoutslll996f ) to detect 
and measure the photometry of sources in each of our 
reduced images. Low signal-to- noise (S/N) borders were 
masked interactively and detections here were excluded. 
We used a DEBLEND_MINCONT value of 0.003, DE- 
TECT_MINAREA set to 8 pixels, and a convolution with 
a 3-pixel Gaussian filter. We experimented with the de- 
tection threshold, finding that a value of 1.2(t detected all 
potential sources in the data. Deblending is an obvious 
concern in close pair studies. We inspected many de- 
blended sources, comparing them to their counterparts 
in the HST z-band imaging, and verified that the de- 
blending algorithm successfully identifies what appears 
to be separate galaxies and not sub-components of in- 
dividual systems. This problem is less extreme in our 
0'.'5 seeing MOIRCS imaging — which smooths over fine 
sub-structure — than in the higher resolution HST data. 
Additionally, we mitigate the potential problem of over- 
deblending by setting a minimum to the pair separation 
(''sep > ^h~^ kpc) in the analysis that follows, although 
this has little effect on our results. 

To investigate the completeness and photometric un- 
certainty of our data, we inserted numerous fake sources 
of varying magnitudes into the reduced images and then 
recovered them using the same SExtractor parameters 
that were applied to real sources. Fake sources were 
given Gaussian profiles with FWHM values similar to 
point sources. We define the image depth by the magni- 
tude corresponding to a recovery rate of 80%. This depth 
implies a slightly lower detection rate (6 0-70%) for real 
galax ies, which are more extended (e.g.. IConselice et"al] 
I2007D . The scatter in recovered magnitudes gives the 
photometric uncertainty as a function of magnitude. 

2.2. ISAAC Near-IR Imaging 

In GOODS-S, JHKg imaging was obtained by an 
ESO/GOODS projecfio using ISAAC on the VLT at the 

10 Program ID: LP168.A-0485 



ESO Paranal Observatory. The survey details and data 
reduction are described in Retzlaff (in preparation). The 
final version 2.0 reduced and calibrated images are pub- 
licly available on the ESO website and, for the Xs-band 
are similar in quality to our MOIRCS data in GOODS-N. 
Covering GOODS-S required 26 pointings given ISAAC'S 
smaller field of view (2.5' on a side). The Ks-h&ud cov- 
erage amounts to 160 arcmin^ with ~95% overlap of the 
GOODS-S ACS imaging. Images were also obtained in 
the J and H bands. After excluding low S/N image bor- 
ders, small gaps (< 1") separate a few of the ISAAC 
pointings. This is not ideal for companion searches, but 
has a negligible impact on our results. The i^T^-band 
depth varies from 24.1 to 25.2 (AB) with seeing less than 
0'.'6 on all images and typically more like 0'.'5. 

As with the MOIRCS data, we used SExtractor to 
perform photometry on the ISAAC data in all three fil- 
ters. We used a DEBLEND_MINCONT value of 0.0005, 
DETECT_MINAREA set to 8 pixels, and a convolution 
with a 3-pixel Gaussian filter. A combined source list 
was assembled with duplicates flagged in the same way 
as above. The image depth and photometric uncertainty 
was estimated using the same procedure of inserting fake 
sources. We note that these data were also used by 
iRawat et all (|2008f ) in a pair fraction study. 

2.3. Matched Catalogs 

We constructed our i^s-selected sample in GOODS-N 
by combining the detections in all of the MOIRCS images 
into a single source list. Because we designed the field 
layout so that the pointings overlap, duplicate sources 
must be removed from this list. We did this by searching 
for those sources that were within 1" of another source 
but were detected on a different MOIRCS image. The 
duplicate with the larger number of SExtractor warning 
flags (bad pixels, saturation, etc.) or worse photometric 
uncertainty was flagged and discarded from the analy- 
sis. The fuU GOODS-N Xs-band catalog contains 8112 
unique sources over 164 arcmin^. 

We matched this i^s-band source list to a number 
of publicly available catalogs in GOODS-N. Beginning 
with the HST data, we used the public ly available vl.l 
GOO DS-N z-selected ACS catalog" (jGiavaUsco et all 
|2004( ). The search tolerance was 0'.'5 for > 22.8 
and 0'.'7 for Kg < 22.8 to account for possible cen- 
troiding problems for bright sources. Where the ACS 
and MOIRCS data fully overlap, «:!86% of the Kg-hand 
sources have an ACS counterpart. For sources with 
Kg < 22.8, this number increases to «98%. 

We then cross referenced our ifg-band sample to sev- 
eral spectroscopic redshift catalogs. For the Team Keck 
Treasury Redshift Survey^^ fTKRS. Iwi^th et al.l [2004). 
we adjusted the TKRS astrometry by -1-0.37" in declina- 
tion and used a search tolerance of 0"6 for Kg > 22.8 and 
0"7 for Kg < 22.8. Of the 4364 TKRS matches we use 
the 1485 with secure redshifts designated by a "zquality" 
code of three or four. We use similar search criteria to 
match our source list to the publicly available compila- 
tion of spectroscopic redshifts^'^ initiated by the survey 
work in lCowie et all (|2004f) . as well as the 1.4 < z < 3.0 

11 http:/ /archive. stsci.edu/prepds/goods/ 

http:/ /tkserver. keck. hawaii.edu/tksurvey/ 
1'^ http://www.ifa.hawaii.edu/ cowie/hhdf/acs.html 
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spectroscopic survey undertaken bv lReddv et alj (|2006D . 
We also cross referenced o ur catalog with t he spectro- 
scopic redshifts obtained bv lTreu et alJ (|2005f ). Including 
TKRS, these comparisons provide a total of 2109 secure 
spectroscopic redshifts for the GOODS-N sample. 

Finally, to obtain U -band observations in order to im- 
prove photo-z estimates, we cross referenced our catalog 
with the deep ground-based photometry obtained by the 
Hawaii Hubble D eep Field North Survey^^ (HHDF-N, 
ICapak et al.ll20d ^. We downloaded the i?-selected cat- 
alogs and adjusted the HHDF-N astrometry by -O'.'l in 
RA and +(f.'2 in Dec to match the ACS-based astrome- 
try used here. The search criteria were the same as for 
TKRS. 

In GOODS-S, a source catalog based on the ISAAC 
data was constructed with the same methods used to 
build the MOIRCS catalog in GOODS-N. Matching to 
the vl.l GOODS-S z-selected ACS catalog was also car- 
ried out in a similar fashion. The ifj-selected ISAAC 
catalog was correlated to source catalogs obtained from 
the ISAAC J and i/-bands {JH photometry was not 
obtained in GOODS-N). The final ISAAC JsT^-selected 
catalog in GOODS-S contains 9043 unique sources over 
160 arcmin^. 

As with the MOIRCS-based catalog in GOODS-N, we 
matched the ISAAC catalog to several publicly avail- 
able surveys. We downloaded Version 1.0 of the VIMOS 
GOODS/ ADP Spectroscopic Survey^^ (jPopesso et al.l 
[200^ and Version 3.0 of the F0RS2 GO ODS/ADP 
Spectroscopic Survey^^ (|Vanzella et al.l [2008I . and refer- 
ences therein), from the ESO/GOODS project website^^. 
Matched photometry, several more spectroscopic redshift 
surveys, and high quality photo -z estimates are com piled 
in the GOODS-MUSIC sample ( jGrazian et al.l[200l . We 
matched the VIMOS, F0RS2, and MUSIC catalogs to 
our ISAAC i^s-selected catalog using a search tolerance 
of 0'.'6. Only ESO redshifts with quality codes of A or B 
or MUSIC codes less than or equal to 1 (signifying high 
confidence in the given redshift) are used in our analy- 
sis. This provides 1683 secure spectroscopic redshifts for 
GOODS-S. 

To summarize, our initial i^s-selected catalog, com- 
plete to Kg « 23.8, consists of 8112 sources in GOODS- 
N and 9043 in GOODS-S for a total of 17155. Of these, 
14998 have ACS counterparts and 3792 have secure spec- 
troscopic redshifts. In §|4]we build a mass-limited galaxy 
sample from these matched catalogs by imposing a mag- 
nitude cut of Kg < 23.57. The resulting sample forms 
the basis of our close pair analysis. Its properties are 
discussed in § 21 

2.4. Supplemental Photometric Redshifts 

Redshifts are a key ingredient in this work. Not only 
are they necessary for examining how the pair fraction 
evolves with time, but they are also used to reject chance 
projections and confirm true physical pairs. Unfortu- 
nately, spectroscopic redshifts are not always available 

^'^ http://www.astro.caltech.edu/ capak/hdf 

Based on observations made with ESO Telescopes at the La 
Silla and Paranal Observatories under Program ID 171.A-3045. 

Observations were carried out using the Very Large Telescope 
at the ESO Paranal Observatory under Program IDs; 170.A-0788, 
074.A-0709, and 275.A-5060. 

http:/ /www. eso.org/science/goods 



and photometric redshifts are needed in both GOODS 
fields. 

Previous photo-z estimates have been made by several 
authors, though not typically for ivTs-selected samples. In 
GOODS-N, high quality estimates are available from the 
HHDF-N Survey, but these use ground-based multiwave- 
length images that are degraded to the worst seeing — 1'.'3 
in the U -band — in order to improve matched photometry 
which is performed in 3" diameter apertures. We found 
that faint companions with separations r < 3" may not 
be detected or may have cont aminated photome try in 
images with such poor seeing (|Capak et al.ll2004f ). and 
so we do not use the HHDF-N photo-z estimates in this 
paper. 

Instead, for the 4703 sources (out of 6812 within ACS 
coverage) in GOODS-N without spectroscopic redshifts, 
we use the Bayesian Photometric Redshift code (BPZ) 
described in iBenfte j (|2000[ ). In GOODS-N, we include 
the [/-band photometry when available from the HHDF- 
N survey combined with 2" diameter photometry in the 
ACS BViz and MOIRCS bands. Comparing to the 
spectroscopic redshifts available, photo-z outliers (de- 
fined by |zspcc — 2phot| > 1) account for 4% of the BPZ es- 
timates, with a-^z/{i+z) ~ 0.1 once outliers are excluded. 

In GOODS-S, photo-zs are required for 6636 sources 
out of the total 8319 sources within the ACS re- 
gion. Precise photo-zs for 5196 were estimated using 
an ACS z-selected sample combined with careful PSF 
fitting and were released in the GOODS-MUSIC com- 
pilation. ^^xchidingoutliers, which account for a few 
percent. ICrazian et al.l (2006) obtain a photo-z precision 
of (Taz/(i+z) ~ 0.03. We use these high-quality photo-zs 
whenever possible in GOODS-S, but still require an addi- 
tional 1440 photo-zs to provide redshifts for full GOODS- 
S catalog. We again use BPZ with ?7-band photometry 
from the MUSIC compilation combined with the ACS 
and JHKs ISAAC photometry measured in 2" diame- 
ters. Comparing to spectroscopic redshifts in GOODS- 
S, we find an outlier fraction of the BPZ photo-zs of 
^10% and crAz/{i+z) ~ 0.11 with outliers excluded. This 
poorer performance of BPZ in GOODS-S compared to 
GOODS-N — despite the addition of the J and H bands — 
seems to result from the shallower ?7-band photometry 
in GOODS-S as well as a discrepancy in the ISAAC filter 
curves 

In both fields comparisons to spectroscopic redshifts 
with z > 1.5 revealed additional catastrophic photo-z 
failures (with incorrect estimates typically assigned to 
zbpz ^ 0.5). Without redshift comparisons at z > 1.5, 
we would have overestimated our photo-z precision. It 
is worth emphasizing that judgments on the quality of 
photo-z estimates from any survey are limited by the 
nature of the spectroscopic comparison set that is avail- 
able. 



3. GALAXY PROPERTIES 

The ISAAC filter transmission curves released on the ISAAC 
instrument website (which give Xab ~ ^Vcga = 1-761) do not 
match the latest calibrations in the final Version 2.0 release. The 
documentation reports that the throughput has been remeasured 
for Version 2.0, giving Kab - ^S'voga = 1-895. Updated ISAAC 
transmission curves — an input to BPZ — have not been made avail- 
able, however. 
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The deep Ks-selected catalogs described above provide 
a powerful dataset for our galaxy pair analysis. Before 
describing the mass-limited sample we draw from this 
dataset and the methodology we use to study pairs, we 
present our techniques for estimating three key physical 
galaxy properties that allow us to characterize the role 
of merging among different galaxy types. 

3.1. Stellar Masses 

The first of these properties is the stellar mass (M,), 
which is estimated using the Bayesian code described in 
iBundv et afl (|2006D . SED fitting is performed using the 
2" diameter BVizKg photometry and the best available 
redshift. The observed SED of each galaxy is compared 
to a g rid of 13440 models from the iBruzual fc Charlotl 
(|2003f ) population synthesis code that span a range of 
metallicities, star formation histories (parametrized as 
exponentials), ages, and dust content. No bursts are in- 
cluded in our models and only those models with ages 
(roughly) less than the cosmic age at each redshift are 
considered. We use a Chabrier IMF (Chabrier 2003) and 
assume a Hubble constant of 70 km Mpc~^. 

At each grid point, the i^s-band M^,/Lk ratios, in- 
ferred M*, and probability that the model matches the 
observed SED is stored. This probability is marginalized 
over the grid, giving the stellar mass probability distri- 
bution, the median of which is taken as the final estimate 
of M^. The width of the distribution provides the un- 
certainty which is typically 0.1-0.2 dex. This is added in 
quadrature to the Ks-hanA magnitude uncertainty to de- 
termine the final error on Af * . Stellar mass estimates for 
galaxies with photo-zs also suffer from the uncertainty in 
luminosity distance introduced by the photo- 2: error and 
the poss ibility of catastrophically wrong redshift infor - 
mation ()Bundv et all 120051 : iKitz 'bichler fc Whiti l2007f ) . 
We will take both into account in the analysis that fol- 
lows. 

More broadly, any stellar mass estimate suffers poten- 
tial systematic errors from uncertainties inherent in stel- 
lar population modeling and various required assump- 
tions, such as the form of the IMF. Several papers 
have stressed the importance of treating thermally pul- 
sating asymptotic gi ant branch stars (TP-AGBs, e.g., 
Maraston et al.ll2006l) an element that is missing in the 
Bruzual & GharlotJ (|200"3 h models. The recent and thor- 
ough investig a tion o f population synthesis modeling by 
[Conrov et al.l (l2008f). however, argues th at M/L ratios 
estimated from IBruzual fc Charlotl (|2003D are largely re- 
sistant to the uncertain contribution from TP-AGBs as 
well as other limitations of current models. Still, it is im- 
portant to recognize that estimates may be affected 
by unrecognized systematic uncertainties at the 0.1-0.2 
dex level. 

3.2. Rest-Frame Color 

We use the inferred rest-frame color to roughly split our 
sample in to blue, star-forming galaxies and red, quiescent 
ones (see IWillmer et al.l 120061 ) . We use the ACS rest- 
frame {B — i) color, with /c-corrections determined from 
the SED fit performed by our stellar mass estimator. For 
each galaxy, the color of every model galaxy in the grid 
is weighted by the probability that the model matches 
the observed SED. The weighted colors are then summed 
over the grid to derive the probability distribution as a 



function of {B — i) color. As with M* we take the median 
of the distribution as the final color estimate. 

The {B — i) distribution across the sample is bimodal, 
as expected. We define red galaxies to be those red- 
der than B — i = 1.3 (AB units) in this filter sys- 
tem. This threshold results in very similar red/blue 
distributions and nu mbe r densities as those reported in 
IWillmer et aP (2006") and IBundv erall (|2006h using data 
from the DEEP2 Redshift Survey. We note however that 
while a restframe color partition is easy to implement, it 
is still a blunt method for discriminating star-forming 
from passive galaxies. For example, contamination from 
dusty star- f orming galaxies o n the red sequence is at least 
-10% fsee lYan et all I2006D and a comparison to mor- 
phological types (see below) for a subset of galaxies with 
0.6 < z < 0.8 and Mq > 10^° shows that ~12% of blue 
galaxies are classified as spheroidal while ~24% of red 
galaxies are classified as disk-like. 

3.3. Morphology 

Galaxy morphology is the final property we consider 
and is important because morphological evolution may 
be driven in part by merging. Morphological clas sifica- 
tion in the GOODS fields was presented in ,Bundv et al.l 
()2005f ) and is publicly available^^. The classification was 
carried out visually to a magnitude limit of zab = 22.5 
which increases the effective mass limit of subsamples 
split by morphology, as we discuss below. Each source 
was assigned one of the following morphological types: - 
2=star, -l=compact, 0==E, 1=E/S0, 2=S0, 3=Sab, 4=S, 
5=Sc d, 6=Irr, 7=Un class, 8=Merger, 9=Fault. As 
m Bu ndv et al.l (|2005f ). we group these definitions into 
the broader categories of spheroidals (types 0-2), disks 
(types 3-5), and irregulars (types 6-8). Visual morpholo- 
gies are available for 1369 sources in GOODS-N and 1200 
in GOODS-S. 

4. METHODS FOR COUNTING GALAXY PAIRS 

In any study of galaxy mergers, it is important to 
distinguish major from minor mergers, since major 
mergers — typically defined to have mass ratios between 
1:4 and 1:1 — have the potential to ra dically affect its 
morphology and star formation (e.g., iNaab fc BurkertI 
l2003h . With many previous studies lacking near-IR data, 
however, pair samples comprised an unknown mix of 
mass ratios. In fact, pairs defined using an optical magni- 
tude differ ence may lead to in flated pair counts from trig- 
gered SF (jBundv et al.ll2004D . an effect that may worsen 
at higher redshifts where the global SFR increases. Thus 
an important goal of this work is to use our i^j-selected 
sample to select only those pairs with stellar mass ratios 
consistent with major mergers. 

We implement this requirement by selecting pairs with 
a magnitude difference^" of AKg < 1.5. Using pairs 
where mass estimates are available for both members, 
we find that this magnitude criterion results in mass 
ratios Mcompanion/Mhost > 1/4 for -80% of the pairs. 
This number is greater than 1/6 for 90%. Pairs are 

http://www.astro.caltech.edu/GOODS_morphs 
Recognizing the assumptions inherent in estimating M«, our 
decision to use the observed magnitude difference as opposed to the 
inferred mass difference is designed to facihtate future comparisons 
to our results. 
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only included if found within a projected annulus of 
5h^^ kpc < Tsop < 20/i~^ kpc determined by the redshift 
of the primary or "host" galaxy, given our assumed cos- 
mology (i?o=70 km s-i Mpc-\ riA/=0.3 and f)A=0.7). 
We avoid double counting by identifying only those com- 
panions that are fainter than their host galaxy. In or- 
der to infer the true physical pair fraction, we use two 
methods to remove chance projections from the raw pair 
counts defined above. The first uses a correction based 
solely on the sky density of background sources while the 
second employs redshift information to identify pairs in 
3D space. 

Before describing these methods in detail, we begin 
with a description of our mass-limited sample of host 
galaxies. It is important that the sample be clearly de- 
fined because the definition of the pair fraction, /pain is 
determined with respect to such a sample. With our defi- 
nition, the number of paired systems is A'pair — /pair-^gai, 
where A'gai is the number of host galaxies. Note that 
this definition differs from that of other authors in that 
/pair is not the fraction of galaxies located in pairs (often 
termed Nc) but instead reflects the potential number of 
mergers within a sample (/pair ~ Nc/2). Finally, because 
our sample allows us to determine not only the fraction, 
but the comoving number density of close pairs, we also 
describe our technique for estimating that quantity in 
this section. Converting the observed pair fractions into 
merger rates requires additional and uncertain assump- 
tions, and we save this topic for § El 

4.1. The Host and Companion Galaxy Sample 

One of the primary advantages of our dataset is that 
both host and companion galaxies are drawn from a well- 
understood, i^s-selected catalog that can be character- 
ized in terms of mass completeness. Based on the differ- 
ent depths of our near-IR images, we optimize our sample 
by choosing a i^s-band magnitude limit of A's = 23.57. 
All companions and host galaxies must be brighter than 
this limit. Because companions are defined to be no more 
than 1.5 magnitudes fainter in the Ks-hand than their 
corresponding host galaxy, this effectively means that our 
host sample is limited to Kg < 22.07. 

Throughout this paper we will adopt three redshift in- 
tervals: 0.4 < z < 0.7, 0.7 < z < 0.9, and 0.9 < z < 1.4. 
Given the Kg-hand limit imposed, the host sample is 
complete in each redshift bin for masses greater than 
logM*/MQ = 9.6, 10.1, and 10.4. The corresponding 
completeness limits for potential companion galaxies at 
these redshifts are log M^/Mg = 9.0, 9.5, and 9.8. When 
morphological samples are considered, an additional z- 
band magnitude limit of zab < 22.5 is effectively im- 
posed since reliable morphologies are difficult to obtain 
at fainter magnitudes. This increases the mass complete- 
ness limits of host galaxies to log Mh,/M0 = 10.3, 10.6, 
and 11.2. Below these values, some hosts may not be 
included in the sample. We highlight the impact of these 
limits on our results as we present them. These stellar 
mass limits are estimated by considering the detection 
efficiency of redshifted SED templates of burst-like stel- 
lar populations with a formation redshift of Zform = 5- 
Because of their high M,/L ratios, these models provide 
conservative and accurate mass limit estimates. 

Since companions are defined by being near their host 
galaxy, the field geometry of the sample is also important. 



We define contiguous survey regions within both GOODS 
fields where the Kg-band detection limit is greater than 
23.57 and the ACS tiles maximally overlap (we exclude 
the triangular, low-S/N perimeters of the GOODS foot- 
print). All of the ISAAC images in GOODS-S meet this 
requirement, giving a search area of 147 arcmin^. We 
exclude three of the 16 MOIRCS images in GOODS-N 
because they are shallower than Kg — 23.57. This gives 
a final area of 139 arcmin^ in the northern field, or 286 
arcmin^ with both fields combined. Only host galaxies 
separated by a distance greater than rmax from the bor- 
der of the defined region are included (this removes ~3% 
from the sample). Sta rs, as deterin ii ied vi sually in the 
morphology catalog of iBundv et al.l (|2005D and also as 
identified spectroscopically, are removed. The final host 
sample includes 2994 primary galaxies, more than 60% 
of which have spectroscopic redshifts. 

Unfortunately, spectroscopic redshifts are not always 
available for both the host galaxy and potential compan- 
ions. We must therefore make statistical corrections that 
account for the contamination of the true physical pair 
fraction from close, line-of-sight projections. We describe 
the two approaches we adopt for this correction below. 

4.2. Method I: Projected Field Correction 

Our first method for correcting the raw pair counts 
is common in pas t studies at high redshift (e.g., 
iLe Fevre et"ani2000l ) and assumes no redshift informa- 
tion for fainter companions is available. We require only 
that potential companions be detected in the X^-band 
and that they satisfy the criteria described above. The 
field contamination for each host galaxy is measured 
according to the i^^-band source density (i.e., number 
counts) observed in the defined survey region in the mag- 
nitude interval fsThost < i^companion < ^i^host + Aif, with 
AK =1.5. Because of cosmic variance, it is important to 
make this measurement separately for the host sample in 
GOODS-N and GOODS-S. For a given host galaxy, the 
measured sky density of possible contaminants is mul- 
tiplied by the area of the corresponding search annulus 
to determine the contamination rate per galaxy. This is 
then summed across the potential host sample to esti- 
mate the total number of field contaminants, A^corr- We 
estimate the number of true physical pairs by subtracting 
A'corr from the total number of potential companions ful- 
filling the search criteria described above. As discussed 
in detail in ^ this method identifies 242 pairs within 
our sample, ~160 of which are expected to by chance 
projections. Uncertainties are determined by adding the 
Poisson error from TVcorr and from the total raw pair 
counts in quadrature. We expect an additional contribu- 
tion from cosmic variance in the relatively small GOODS 
fields. 

4.3. Method II: Probable Redshift Confirmation 

The vast majority of faint companions identified in 
the /fs-band catalog are also detected in the ACS and 
ground-based optical imaging which means that photo- 
metric redshifts can be estimated for these sources. A 
sizable fraction of companions (~30%) also have spec- 
troscopic redshifts. In our second pair counting method, 
we use the redshift information available for faint near-by 
neighbors to find those that are likely to be physically as- 
sociated with their brighter host. This amounts to an ad- 
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ditional redshift constraint on the identification of com- 
panions. For both the potential hosts and companions 
we define a redshift uncertainty of cr^ — 0.002 for galaxies 
with spectroscopic redshifts, cr^ — 0.08(1 + z) for galax- 
ies with photometric redshifts^^, and an unbounded cr^ 
for the remaining few galaxies without optical detections 
and therefore with no redshift information. When the 
redshift difference satisfies Az^ < a^ f^^^^ + al^^^p^^^^^, 
the pair is selected^^ with this method. 

Because of the large redshift uncertainty, some of the 
host-companion pairs defined by photometric redshifts 
will be close in projection but not physically associated. 
As in our first pair counting method, a contamination 
correction is needed although it will of course be smaller 
since most contaminants a re already exclud e d by their 
discrepant redshifts. As in iKartaltepe et all (|2007D . we 
determine this correction factor by randomizing the co- 
ordinate positions of all the sources in our catalog (but 
retaining all redshift and magnitude information). Any 
companions identified in the randomized data are clearly 
chance superpositions. Thus, by repeating this exercise 
for 100 randomized datasets, we can estimate the average 
contamination rate, its uncertainty, and the dependence 
of these quantities on A/* and pair type (see below). Note 
that in the limit that no redshift information is available 
for faint neighbors this approach is equivalent to the first 
method described above. However, for surveys that con- 
tain dense environments like galaxy clusters (GOODS is 
too small to sample such environments), this technique 
underestimates the needed correction because it does not 
account for the increased rate of chance projections in 
dense regions. 

Finally, we note that the possibility of catastrophic er- 
rors in photo-z estimates will cause some true physical 
companions to be disregarded as chance projections. We 
can determine the importance of this effect by examin- 
ing the quality of our various photo-z estimates (§ 12. 4p 
and estimating the likely number of catastrophic photo-z 
failures that cause us to underestimate the true number 
of pairs in method II. Based on this analysis we correct 
a-U /pair values derived from method II upwards by 1%. 
Using method II, we find 89 pairs, ~ 41 of which are 
expected to be chance projections. 

4.4. Estimating Comoving Number Densities 

One of the important quantities we wish to determine 
in this work is the comoving volumetric merger rate as a 
function of the stellar mass of the host galaxy. As we dis- 
cuss in § [6l this "merger rate mass function" provides a 
powerful tool for understanding the role of galaxy merg- 
ers in driving mass assembly as well as evolution in the 
numbers of different populations such as ellipticals and 
red-sequence galaxies. 

We estimate such ma.ss fun ctions using the simple 
Vmax technique ()Schmidtj|1968D . We weight host galaxies 
and associated pairs by the maximum volume in which 
they would be detected within the Xs-band limit in a 

The value of 0.08(1 + z) is somewhat arbitrary and represents 
a compromise between the BPZ photo-2 uncertainty {~0.1) and 
the MUSIC photo-2 uncertainty (~0.03). 

^■^ For the pairs in our sample identified by method II, ~26% 
feature spectroscopic redshifts for both members, ~39% require at 
least one photo-2 ~28% require photo-zs for both members, and 
~7% are without redshift information for the fainter companion 



given redshift interval. In practice, our magnitude lim- 
its impose restrictions only in the highest redshift bin 
(0.9 < z < 1.4). For each host galaxy i in the redshift 
interval j, the value of V^^^ is given by the minimum 
redshift at which the galaxy would drop out of the sam- 
ple, 

^Ca. = dn—dz (1) 

where dQ is the solid angle subtended by the survey area, 
and dV/dz is the comoving volume element. The redshift 
limits are given as, 

Zhigh = min(z^„^ '^KuJ 

Zlow = ^miri (3) 

where the redshift interval, j, is defined by [z;^j„, z^j^j^.] 
and z , refers to the redshift at which the galaxy would 
still be detected below the Ks-hand limit for that partic- 
ular redshift interval. We use the best-fit SED template 
as determined by the stellar mass estimator to calculate 
z]^^ , thereby accounting for the /c-corrections necessary 
to compute Vmax values (no evolutionary correction is 
applied) . 

5. RESULTS 

5.1. The Stellar Mass Dependent Pair Fraction 

As discussed in § [TJ the hierarchical framework that 
underlies current galaxy formation models argues that 
the mass scale on which galaxies assemble increases with 
time. Observations of star formation quenching and mor- 
phological evolution show that these processes are set 
by a mass scale that decreases with time (downsizing). 
Meanwhile, galaxy mergers are thought to be key drivers 
of both hierarchical assembly and the evolution of SF 
and galaxy morphology. So, understanding the mass de- 
pendence of merging provides a critical test of the role 
mergers play in galaxy mass assembly and evolution. 

We begin our investigation of how galaxy mergers de- 
pend on mass by presenting the pair fraction measured 
using the two methods described in the previous sec- 
tion. With method I, in which no redshift information 
for companion galaxies is used, 242 pairs are found, ~160 
of which are expected to be chance projections. The 
method I pair fractions, partitioned by mass and redshift, 
are plotted in the top panel of Figure[TJ With method II, 
which takes companion redshifts into account, 89 pairs 
are found, ^41 of which are expected projections. The 
resulting method II pair fractions are plotted in the bot- 
tom panel of [1] In both cases, pair fractions are plotted 
for three redshift intervals: 0.4 < z < 0.7, 0.7 < z < 0.9, 
and 0.9 < z < 1.4. The full sample (X symbols) as well 
as the dependence on the mass of the host, as indicated 
by the differently sized circular symbols, is shown. All of 
the plotted points are complete in terms of stellar mass 
except for the highest redshift 10 < logM*/M0< 10.5 
bin which is ~80% complete as a result of the ^g-band 
limit. The results from both methods are listed in Table 

m 

We first note that the two panels in Figure [1] are sim- 
ilar, reinforcing the utility of the two methods. Both 
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Fig. 1. — Contamination-corrected fraction of paired systems 
with 5h~^ kpc < r^cp < 20h~^ kpc among galaxies in three red- 
shift intervals. The top panel shows results determined using 
method I, the bottom panel method II. Abscissa values have been 
spaced to keep the data points from overlapping. The X symbols 
denote the full sample of primary (host) galaxies while the size of 
the circular symbols indicates the result of dividing these into stel- 
lar mass bins, as labeled. Predic ted v alues of /pair ba sed o n the 
2-point correlation function (see § 16.11 1 of Zeh avi et al. (2005!) (us- 
ing SDSS) are overplotted in the low-z bin of the to p panel while 
the mid-z bin shows predictions based on Coil et at| | |200^ ( us- 
ing DEEP2). The mass dependence of the clustering predictions 
is based roughly on the luminosity range probed by the clustering 
analyses. 



plots show that roughly 2-5% of galaxies in a mass- 
limited sample host fainter companions. The value 
of /pair depends on the maximum allowed pair sepa- 
ration, Tniax- The average value of /pair increases by 
an amount of 0.03 for r^ax = 25h^^ kpc compared to 
rmnK = 15/i~^ kpc. This difference does not impact the 
derived merger rate because the exp ected merger effi- 
cienc y decreases with larger r^ax fe.g.. lPatton fc Atfieldl 
I2008D . We use 5/i"^ kpc < < 20/1"^ kpc in what 
follows. 

The key result in Figure [1] is the increase in /pair among 
more massiv e host ga l axies. Hints of this trend were 
observed by iXu et all (j2004) in their iiT-selected study 
of close pairs at {z) = 0.03. We confirm it here with 
moderate significance to 2; > 1 using both pair count 
correction methods. We verified that these results are 
still apparent for different values of rgep including rgop = 
15/i~^kpc and Vsc-p = 25/i~^kpc. The mass dependence is 
also apparent in each GOODS field separately, suggesting 
that this trend is not a result of cosmic variance. 

As we discuss in § [51 some of the increase at higher 
masses may reflect the stronger clustering of such galax- 
ies. In that section, we return to the important ques- 
tion regarding the extent to which the mass dependence 
observed for the pair fraction translates into a mass de- 
pendence for the merger rate. Still, the fact that the 
lowest mass bins in Figure [T] are nearly consistent with 
a zero pair fraction implies that the merger rate inferred 
from this and other studies is dominated by higher mass 
galaxies. We show in § 16.51 that this leads to the result 



that major mergers are unlikely to be the sole mechanism 
behind the formation of spheroidal galaxies and the red 
sequence. At the same time, the higher merger rates im- 
plied for massive galaxies shows that they are continuing 
to assemble after z \. This is an important piece of 
direct evidence for hierarchical growth which we seek to 
quantify in § 16.51 

Finally, we comment on the possibility that /pair 
evolves with redshift. Averaging over the full mass range, 
Mf, > IO^^Mq, we find no statistically significant evolu- 
tion. With /pair oc (1 -I- z)™, we measure m = 1.6 ± 1.6 
using method I and m — 0.3 ± 1.4 using method II. 

5.2. The Dependence on Host Galaxy Color and 
Morphology 

In addition to their role in mass assembly and their 
potential to affect the morphology and SFR of galaxies, 
mergers have also been invoked as a means to build the 
massive end of the red sequence. Successive dry merg- 
ers between early-type galaxies at lower n iasses can in- 
crease the numbers at h igher masses (e.g.. Ivan Dokkuinl 
l2005t iFaber eTall 12001 . In this scenario, mergers must 
clearly be common within the red-sequence. We can be- 
gin to test this by characterizing the dependence of /pair 
on host galaxy type. Our goal is not only to count the 
frequency of dry mergers, but also the remaining fraction 
of wet mergers since only these events can be responsi- 
ble for morphological transformation or star formation 
quenching. 

Both pair counting methods offer insight on the ques- 
tion of how the pair fraction varies with galaxy type. 
Since the first method relies on an entirely statistical 
field correction, we can only use it to study trends with 
host galaxy type. The properties of the companion and 
host together — the wet and dry merger frequency, for 
example — are accessible with method II and discussed 
below. Method I allows us to divide the host sample by 
both mass and type. Figure [2] shows the result when host 
galaxy type is defined using the color bimodality to sep- 
arate star-forming from passive systems (also see Table 
[1]) . Bins corresponding to the most massive blue galaxies 
have not been plotted because they contain fewer than 
10 hosts. 

While the uncertainties have increased compared to 
FigurelU there is weak evidence that the mass-dependent 
trend observed in Figure [T] is reflected in the star- forming 
properties of the host galaxies, with higher pair frac- 
tions found for quenched, red hosts, in large part because 
such galaxies dominate at the highest masses. This sug- 
gests a role for dry mergers in building the most mas- 
sive early-type galaxies. However, a second interpreta- 
tion discussed further in § [6] is that Figure [2] may re- 
veal the greater degree of clustering among massive red 
galaxies. Finally, we note that the number of red hosts 
with companions in the lowest mass bin is consistent with 
zero. If confirmed with future studies, this observation 
would indicate that major mergers among red systems 
with < IO^^Mq cannot contribute significantly to the 
increasing abundance of red sequence galaxies at higher 
masses. Transforming blue galaxies into red objects via 
quenched star formation would therefore appear to be a 
more important mechanism for red sequence growth at 
< 10"Mo. 

In Figure [3] we explore the relationship between pair 
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Fig. 2. — Color and mass dependence of the field-corrected pair 
fraction, /pair using method I. The host sample has been divided 
into "blue" (blue triangles) and "red" (red squares) based on rest- 
frame color. Symbol size indicates stellar mass as in Figure [l] 
Bins corresponding to the most massive blue galaxies have not 
been plotted because they contain fewer than 10 hosts. 

fraction and galaxy morphology. Here the host sample 
has been restricted to zab < 22.5 which makes some 
mass bins incomplete (although we note we are mea- 
suring fractions not absolute quantities). Bins that are 
less than ~80% complete incomplete are indicated by 
lighter colors and open symbols and will tend to have 
an artificially lower number of redder galaxies. Given 
the large error bars, any trends observed here will re- 
quire confirmation from larger studies. Still, there is 
a hint that spheroidals and irregulars host more com- 
panions than disk galaxies. Spheroidals, which are ex- 
pected to be mostly red and quiescent, also tend to 
exhibit higher pair fractions among high-mass galaxies 
with logM,j/M0> 10.5 (although this is not the case at 
z - 0.5). 

One intriguing and more robust result of Figure 
[3] is the appearance of a high fraction of irregular 
hosts at intermediate masses. In some ways, this 
is unexpected because irregulars represent only 10- 
20% of the galaxy popula tion in our redshift range 
with logMJMfT, > 10 (Brinchmann fc EUid l2000l : 
iBundv et al] [200a [Pmllelia et al. 2006). They should 
therefore be less likely to host companions, if pairs were 
randomly distributed among all galaxy types. Even in in- 
complete bins which are biased against red spheroidals, 
host galaxies are more likely to have irregular morpholo- 
gies than disk-like ones. Because we ignore compan- 
ions with Tsop < 5/i~^kpc, the prevalence of irregu- 
lar hosts tentatively suggests that even more distant 
companions can imprint significant morphological dis- 
turbances on their host galaxy, perhaps after executing 
a first pass (e.g.,OLotz ct al. 2008a). A similar conclu- 
sion was reached by Li et al. (2008) using SDSS. They 
also find an enhancement of SF among close pairs which 
could lead to morphological pecuharities. It stands to 
reason that not all irregular systems are in the last stages 
of a significant merger, which may in part explain why 
merger rates based on disturbed morphologies are typ- 
ically higher than those inferred from pair counts (e.g., 
lConselicel[2008l : iLotz et al.ll2008H ). 

Taken together, the type dependent pair fractions we 
have presented suggest that galaxies with higher masses, 
redder colors, and early- type morphologies tend to host 
companions more frequently than their lower mass, star- 
forming, and disk-like counterparts. This lends support 



Fig. 3. — Morphological dependence of the field-corrected pair 
fraction using method I. The host sample has been restricted to 
^AB < 22.5 and visually classified by ACS morphology into disks 
(blue triangles) and spheroidals (red squares) and irregulars (ma- 
genta stars). Symbol size indicates stellar mass as in Figure[T]with 
the smallest samples denoting the full mass range. Data points 
corresponding to bins that are significantly incomplete in stellar 
mass are indicated by lighter colors and open symbols. 

to the concept of dry mergers as a mechanism for red se- 
quence growth, especially at the highest masses, a point 
we return to in § 16.41 However, dry mergers do not ap- 
pear to be important for building the red sequence at 
M» 10^°Mq because the number of red hosts at these 
masses is nearly zero. At the same time, the enhanced 
pair fraction for more massive red and spheroidal galax- 
ies may simply reflect the denser environments of these 
galaxies. The question of whether the derived merger 
rates are also enhanced in such environments is discussed 
in §[6l 

5.3. The Frequency of "Dry" and "Wet" Pairs 

In the top panel of Figure[4]we use method II combined 
with the rest-frame color of pair members to distinguish 
the fraction of pairs expected to result in dry and wet 
mergers (also see Table [2|) . Red-red pairs are identified 
as dry merger candidates. Blue-blue and mixed pairs 
indicate the presence of cold gas reservoirs and therefore 
represent potential wet mergers. For the few secondary 
(companion) galaxies without a redshift, the pair type 
cannot be determined and is labeled as "N/A." Note that 
a statistical correction as described in § 14.31 has been 
subtracted from each category so that negative values 
are possible in some cases. 

An interesting evolutionary signal appears in Figure 
m First, it appears that dry, red-red pairs are less 
frequent than their wet counterparts (the sum of blue- 
blue and mixed) in the highest redshift bin but be- 
come increasingly important with cosmic time, even- 
tually dominating the low redshift bin. At the same 
time the number of blue-blue pairs decreases to zero. 
These trends are significant at the 1-a level and may 
simply reflect the changing nature of the galaxy popu- 
lation with masses log M^/M© > 10.5 which is increas- 
ingly dominated by red galaxies at late times. If our 
sample size allowed us to consider only high mass galax- 
ies (e.g., log(M*/Mo) > 11) we might expect dry pairs 
among such galaxies to show an even earlier increase 
since th is populatio n evolve s into quenched red galax- 
ies fir st (|Bundv et a l. 2006; Borch et al.ll2006l : iBell et all 
|2007| ). However, it is important to emphasize that the 
dry pair fraction is consistent with what is expected given 
the makeup of galaxies from which the pairs are drawn. 
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Redshift 

Fig. 4. — Fraction of dry, wet, and mixed pairs as defined by 
both color (top panel) and morphology (bottom panel), for host 
galaxies with logAf*/M0 > 10.5. Pairs are confirmed using spec- 
troscopic and photometric redshifts, and a correction for contam- 
ination and catastrophic photo-zs has been applied (method II). 
Pairs including secondaries without a redshift are listed as "N/A." 
In the bottom panel, early-type refers to galaxies with spheroidal 
morphologies (E/SO/Sab) while late-type includes both disk-like 
and morphologically irregular galaxies. An additional magnitude 
cut of ZA.B < 22.5 applies to the bottom panel. 



A similar conclusion was reached by iLin et aP (|2008l ) 
who studied the fractions of blue-blue, red-red, and 
mixed systems among kinematic pairs identified in the 
DEEP 2 survey. They also found a higher number of 
blue-blue pairs at z ~ 1 compared to the present day 
combined with an increasing number of red-red and, to 
a lesser extent, mixed systems. While that work benefits 
from spectroscopic redshifts for all pair members, our 
analysis has the advantage of being mass-limited and 
capable of selecting only major pairs. Furthermore, it 
does not require weighting to account for biases in the 
spectroscopic target selection. Given these very differ- 
ent but complementary methods, it is encouraging to 
see agreement both in the derived type-dependent pair 
fraction (note the different definitions of pair fraction: 
/pair — 0.5iVc) and its evolution with redshift. 

Our data also allow us to examine the frequency of 
pairs defined by the morphological type of both mem- 
bers. This is shown in the bottom panel of Figure|3]where 
it is important to remember that in the highest redshift 
bin, 30-40% of red galaxies with log(M,/M0) > 10.5 
will be missed because of the zab < 22.5 morphology 
cut. In this panel we label disk- like and irregular galax- 
ies as "late-type" and spheroidals as "early-type." We 
see no statistically significant trends, although there is a 
hint that late-late pairs are on par with early-early and 
mixed pairs in the highest redshift bin, but become sub- 
dominant in the mid- and low-z bins. Although larger 
datasets will be needed to confirm this, such an effect 
would be expected as the fraction of late-type galaxies 
declines with time. 

As we discuss further in § [6l Figure [4] has important 
implications on the role of wet mergers in transforming 
blue (disk) galaxies into red (spheroidal) galaxies. Both 



TABLE 2 

Method II: Dry, Wet, and Mixed Mergers 



Sample 


z 


Np 


Ncorr Ngal 


/pair (%) 




Color, lof 


iM^IMq > 10.5 






Blue-Blue 


0.4-0.7 


4 


2.3 


514 


1.3 ± 0.5 


Blue-Blue 


0.7-0.9 


7 


3.3 


380 


2.0 ± 0.8 


Blue-Blue 


0.9-1.4 


24 


9.4 


793 


2.8 ± 0.7 


Red-Red 


0.4-0.7 


7 


0.7 


514 


2.2 ± 0.5 


Red-Red 


0.7-0.9 


5 


0.9 


380 


2.1 ± 0.6 


Red-Red 


0.9-1.4 


8 


1.3 


793 


1.8 ± 0.4 


Mixed 


0.4-0.7 


8 


2.6 


514 


2.0 ± 0.6 


Mixed 


0.7-0.9 


10 


3.5 


380 


2.7 ± 1.0 


Mixed 


0.9-1.4 


10 


5.8 


793 


1.5 ± 0.5 


N/A 


0.4-0.7 





1.1 


514 


0.8 ± 0.2 


N/A 


0.7-0.9 


1 


1.0 


380 


1.0 ± 0.4 


N/A 


0.9-1.4 


1 


1.2 


793 


1.0 ± 0.2 


Morphology, logAf./MQ > 10.5 


Late-Late 


0.4-0.7 


2 


1.5 


482 


1.1 ± 0.4 


Late-Late 


0.7-0.9 





0.4 


261 


0.8 ± 0.3 


Late-Late 


0.9-1.4 


3 


0.5 


322 


1.8 ± 0.6 


Early-Early 


0.4-0.7 


4 


0.4 


482 


1.8 ± 0.4 


Early-Early 


0.7-0.9 


3 


0.2 


261 


2.1 ± 0.7 


Early-Early 


0.9-1.4 


3 


0.0 


322 


1.9 ± 0.5 


Mixed 


0.4-0.7 


7 


1.4 


482 


2.2 ± 0.6 


Mixed 


0.7-0.9 


4 


0.6 


261 


2.3 ± 0.8 


Mixed 


0.9-1.4 


1 


0.3 


322 


1.2 ± 0.4 


N/A 


0.4-0.7 


4 


1.5 


482 


1.5 ± 0.5 


N/A 


0.7-0.9 


1 


1.1 


261 


1.0 ± 0.6 


N/A 


0.9-1.4 


2 


0.9 


322 


1.3 ± 0.5 



Note. — /pair has been increased by 1% to account for 
losses due to catastrophic photo-z errors. 

indicate that the implied number of such wet mergers is 
a factor of ~2 less than the total merger rate. These 
figures also show that the rate of dry merging is not en- 
hanced beyond what is expected given the makeup of 
galaxies in a mass-limited sample. In § [6] we show ex- 
plicitly that neither morphological transformations oc- 
curring after wet mergers nor mass buildup caused by 
dry mergers appears sufficient to explain the increasing 
number of intermediate to very massive early-type galax- 
ies. 

5.4. Comparison to Previous Work 

In this subsection, we discuss our results in the con- 
text of previous work. A number of close pair stud- 
ies have been carried out us in g local surveys with 
z < 0.3 fe.g iXu et all [200l ID c Propris et all 120051: 
d] l2006l : iPatton fc Atfield .2 



iMasiedi et all 



,20q|). These tend 



to find very low pair fractions, often only a few percent. 
Using our full sample with log M^/M© > 10, we find 

/i 



pair 



0.03 ± 0.02 at z - 0.5. This is slightly above 
but consistent wi t h the value of /pair — 0.01 reported by 
iPatton fc Atfieldl (|2008f) who analyzed a volume-limited 
r-selected SDSS sample (note that /pair ~ O.bNc). They 
selected pairs with r-band luminosity ratios greater than 
1/2 which is more stringent than our « 1/4 ratio 
threshold and should lead to lower /pair since not al l ma- 
jor mergers are selected. iPatton fc Atfieldl (|2008f ) find 
no luminosity dependence in their pair fraction estimate 
(with Tsep =5-20 kpc) in contrast to the mass de- 
p endence found h ere . 

IXu et all (j2004[ ) use the combination of the 2DFRGS 
and 2MASS surveys to measure both the fraction and 
number density of close pairs drawn from a X-selected 
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sample at {z) — 0.03. With a magnitude threshold of 
Aif < 1, their study is an excellent low-z analog to ours, 
and they estimate /pair ~0. 02-0. 03 (their pair fraction is 
also defined as 2 /pair) in good agreemen t with our lowest 
redshift bin. As mentioned previously, IXu et all (|2004l ) 
also find weak evidence for a higher value of /pair for more 
massive galaxies. Our analysis confirms this finding and 
demonstrates that it continues to z 1. 

Studies at higher redshifts tend to focus on the evolu- 
tion of the pair fraction with redshift as parametrized by 
TO. This number is poorly cons trained by recent work , 
with TO varying from to 4 (see iKartaltepe et al.ll2007h . 
Leaving evolution aside, better agreement in the average 
merger rate is being achieved, however. For example, 
our new results are consistent wit h the initial near-IR 
st udy undertaken by iBundv et al.l ([2004) and the work 
of lRawat et al.l (|2008[ ) who perform a similar analysis us- 
ing only spectroscopic redshifts in GOODS-S. 

A number of optical studies of the pair fraction also 
deliver similar results, although some caution is needed 
in making comparisons. In addition to the bias in opti- 
cal samples due to SF, our present work reveals a second 
effect. Depending on the selection method, optical sam- 
ples may trace different stellar mass ranges as a function 
of redshift. Since we have shown that the pair fraction 
depends on Af,, this can bias the pair fraction at dif- 
ferent redshifts. With this in mind, we not e that the 
value o f /pair =0.02-0.03 measured by both iBell et all 
(I2006ar) using COMBO-17 as w ell as the k inema tic pair 
study of iLin et al.l ()2004f ) and iLin et all (|2008f ) agrees 
with our findings here. Similar values were also found by 
IKartaltepe et a l. (2007) and Hsi eh et al.l ()2008[ ) and, for 
z < 0.8. bv ldTRavel et al.. (20081 )! 

Not all merger studies find such low values, however. 
Our pair fraction is lower by a factor of 3 -4 and shows 
less evolution compared to iLe Fevre et al.l ([lOOO) , likely 
as a result of the SF bias. At the same time, morpho- 
logical derivations of the merger fraction tend to find 
values of /pair « 0.1 at 0.4 < 2 < 1.4, 2 -3 times higher 
than the results o f our pair analysis fe.g. . lConselicel[2003l : 
iLotz et al.|[2008b[ ). This may reflect the influence of non- 
uniform star formation, minor mergers, and even flybys 
on galaxy mor phology as discus sed with respect to Fig- 
ure [3] (also see lLotz et al.ll2068al) . 

6. DISCUSSION: THE GALAXY MERGER RATE 

6.1. Insight from Galaxy Clustering 

Our goal in measuring the pair fraction is to constrain 
the rate at which galaxies merge, but doing so requires 
knowing the fraction, C,ng, of close pairs that will even- 
tually merge and on what timescale, Ting, they coalesce. 
Both quantities are uncertain, but previous work has typ- 
ically assumed that roughly half of close pa irs will merge 
with a typical timescale of 0.5 Gyr (e.g., iPatton et al.l 
|2000[ ). Under these assumptions, the pair fractions pre- 
sented here would indicate that mergers occur more fre- 
quently among massive, red galaxies as compared to their 
less massive and bluer counterparts. But if either Cmg or 
Tmg depend on mass or galaxy type, this conclusion could 
be incorrect. As we now show, evidence from galaxy clus- 
tering may imply such a dependence. 

The galaxy 2-point correlation function is typically 
measured on scales larger than 100 /i^^ kpc and is usually 
fit as a power law ^ = (r/rg)^^. By extrapolating this fit 



to the radii of interest for the close pairs in this paper, 
one can estimate the pair fraction predicted by clustering 
measurements made on larger scales. Assuming f ^ 1 
and using our definition of /pair, 



,air ~ 47rng 



.(1+^) 



r'^£^{r)dr, 



Ann. 



3-7 



^r^oi^ + z)'~''{rl-Z-r'^,^^), 



(4) 
(5) 



where rig is the comoving galaxy number density of po- 
tential companions, and the (1 -I- z) factor accounts for 
the fact that rmin and rmax are defined in physical and 
not comoving coordinates. 

We can use estimates of ro and 7 from various cluster- 
ing analyses to see how well the clustering signal extrap- 
olated to small scales agrees with the pair fractions we 
measure. For example, IZehavi et al.l ()2005f ) measure lu- 
minosity dependent clustering in the SDSS. Taking their 
three brightest luminosity bins to roughly correspond 
to our three stellar mass bins, we use ro = 5.52, 6.16, 10.0 
in units of /i-i Mpc with 7 = 1.78,1.85,2.04. We 
approximate Uq by integrating the mass functions of 



iBundv et all ()2006l ) and use Ug = 4,4,0.8 in units of 
10~^ Mpc~^ with h = 0.7. With these values and set- 
ting z = 0.5 (note that SDSS has an average z w 0.1), we 
find a predicted /pair of 0.03, 0.06, and 0.10 for the three 
"mass" bins. These are overplotted in the first redshift 
bin of the top panel in Figure [TJ 

Extending this comparison to higher redshifts is diffi- 
cult because the smaller sizes of high-z samples do not 
adequately probe the most massive or luminous ga laxies. 
Still, we can compare to the z ~ 1 results of , Coil et al.l 
(|2006[ ) where we take ro = 3.78/1"^ Mpc and 7 = 1.68 
from their median Mb ~ —20.3 sample 1 to corre- 
spond to our lowest mass bin (log M^/M© =10-10.5) and 
ro = 4.21/i"i Mpc with 7 = 1.9 from their Mb « -21 
sample 4 for our second mass bin (log M^/M© =10.5-11). 
We find /pair «0.01 and 0.04 and these results are plot- 
ted in the middle redshift bin of the top panel of Figure 
[TJ It should be emphasized that Mb is not a good tra cer 
of Af, and that the s ample used bv ICoil et al.l (|2006D is 
not mass-limited (see lCoil et al.ll2008| ). 

While the comparison between /pair and predictions 
from the 2-point correlation function is clearly approx- 
imate, the resulting agreement is striking. Indeed, the 
lower corre lation lengths of blue versus red galaxies (e.g., 
ICoil et aL|[200S ) also seem to be consistent with the lower 
pair fraction of blue hosts shown in Figure [H 

This suggests that the connection between /pair and 
clustering is a strong one and that the trends observed 
here may be expected given the way the correlation func- 



for red, early-type systems (e.g., 
IZehavi et al.ri 2q05t iLi et al.l [200a 



Norbere et al. 


2002 


Meneux et al. 


2006 



Coil et al. 20061 l2008t I Meneux et al.l I2008D . In other 
words, the trends we observe with /pair may be driven 
by regions of higher density which tend to host more 
massive, red-sequence systems. 

6.2. Estimates of the Merger Timescale from 
Simulations 
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TABLE 3 

Fractional Merger Rates, TZn 



Patton . 



Kitzbichler . 



logM* 
(log/ifoMo) 



/pair 

(Method I) 



(Gyr) (Gyr-i) 



7^n 



(/ifo^ Gyr) (h70 Gyr-i) 



0.4 < 


z 


< 0.7 


10.0 < logM* 


0.03 


0.58 


0.5 


0.036 


1.0 


2.0 


0.016 


0.7 < 


z 


< 0.9 


10.0 < logM. 


0.05 


0.58 


0.5 


0.060 


1.0 


2.0 


0.026 


0.9 < 


z 


< 1.4 


10.0 < log Af, 


0.06 


0.58 


0.5 


0.068 


1.0 


2.0 


0.029 


0.4 < 


z 


< 0.7 


10.0 < k 


)gM. < 10.5 


0.03 


0.55 


0.5 


0.034 


1.0 


1.9 


0.016 


0.7 < 


z 


< 0.9 


10.0 < Ic 


)gM* < 10.5 


0.05 


0.55 


0.5 


0.057 


1.0 


1.9 


0.027 


0.9 < 


z 


< 1.4 


10.0 < k 


)gM* < 10.5 


0.03 


0.55 


0.5 


0.038 


1.0 


1.9 


0.018 


0.4 < 


z 


< 0.7 


10.5 < k 


)gM* < 11.0 


0.03 


0.58 


0.5 


0.031 


1.0 


1.4 


0.019 


0.7 < 


z 


< 0.9 


10.5 < k 


)gM* < 11.0 


0.04 


0.58 


0.5 


0.048 


1.0 


1.4 


0.030 


0.9 < 


z 


< 1.4 


10.5 < k 


)gA/, < 11.0 


0.07 


0.58 


0.5 


0.080 


1.0 


1.4 


0.050 


0.4 < 


z 


< 0.7 


11.0 < logAf* 


0.05 


0.46 


0.5 


0.046 


1.0 


1.0 


0.051 


0.7 < 


z 


< 0.9 


11.0 < logM* 


0.07 


0.46 


0.5 


0.067 


1.0 


1.0 


0.075 


0.9 < 


z 


< 1.4 


11.0 < logM* 


0.09 


0.46 


0.5 


0.082 


1.0 


1.0 


0.092 



TABLE 4 

Volumetric Merger Rate Mass Function 





0.4 < 2 < 0.7 




0.7 < z < 0.9 


0.9 < 


z < 1.4 




Stellar Mass Range 


log "I- log If wot 


logo- 


log 'S/ log v^wct log o- 


log 't log 


'fwot 


logo- 


10.0 < logAf. < 10.5 


-3.77 -3.96 


0.31 


-3.69 -3.94 0.25 


-4.16 


-4.37 


0.27 


10.5 < logAf* < 11 


-3.90 -4.09 


0.29 


-3.89 -4.14 0.26 


-3.81 


-4.03 


0.13 


logM* > 11 


-4.51 -4.70 


0.28 


-4.38 -4.63 0.24 


-4.50 


-4.72 


0.17 


Volumetric Formation Rate of Spheroidals 


10.25 < log A/* < 10.55 


-3.18 


0.14 


-3.12 ... 0.09 


-4.50 




0.42 


10.55 < logM. < 10.85 


-3.92 


0.55 


-2.91 ... 0.08 


-3.37 




0.09 


10.85 < logM, < 11.15 


-3.57 


0.33 


-3.02 ... 0.10 


-3.55 




0.16 


log A'/* > 11.15 


-4.06 


0.52 


-3.44 ... 0.17 


-3.60 




0.13 



Note. — The values in this table are plotted in Figure [5] The symbol ^ denotes a rate per unit volume (per 
logarithmic interval) with units of h^Q^Mpc~'^dex~^Gyr~'^. It should be distinguished from <&, often used to 
represent mass or luminosity functions, ^wct indicates the approximate volumetric merger rate after removing 
dry mergers. We use logo- to designate the associated statistical uncertainty (in dex). Stellar masses have units 
2, 



of H^qMq. 

Perhaps the best "way to estimate Cmg and T^g 
is through detailed simulations. Unfortunately, fully 
numerical merger simulations are expensive and can 
only be performed for a handful of systems (e.g., 
iBovlan-Kolchin et al.|[2007[ ). Cosmological simulations, 
mean-while, are limited in resolution and require var- 
ious analytic assumptions regarding dynamical fric- 
tion, tidal interact i ons, and other details. Both 
iKitzbichler fc White! (|200l and lPatton fc Atfieldl (pCiOSh 
use the Millennium Simulation to determine the fre- 
quency and timescales of merg ing pairs, and af- 
ter ac counting for the fact that IKitzbichler fc White! 
(l2008f) absorb Cmg int o their estimate of (Tmerge) (see 
iPatton fc Atfieldl |2008[ ). the two studies derive similar 
r esults on average. 

IKitzbichler fc Whitd (|2008l ) find that (Tmergo) depends 
inversely on stellar mass as M~°-^. This actually en- 
hances the mass dependence of the merger rate beyond 
what we have observed for /pair- In other words, the 



higher /pair for massi ve galaxies transla t es int o an even 
faster merger rate. iPatton fc Atfieldl (|2008f l. on the 
other hand, assume that Tmg, defined for real-space pairs 
destined to merge, is a roughly constant 0.5 Gyr and 
instead investigate the dependence of Cmg (defined as 
/3D in the ir paper) on luminos i ty. As might be expected 
from the .Kitzbichler fc Whitl ()2008D analysis, they find 
increases with luminosity, resulting in more 
efficient merging for more luminous pairs (equivalent 
to a lower value of (Tmorgc) for more massive galaxies). 
However, it is interesting that their most luminous pairs 
{Mr ~ -22.8 or log M*/M0 « 11.2) show a drop in C„g 
which may be evidence of the c l usterin g effect discussed 
above. The IPatton fc Atfieldl (|2008f ) values therefore 
have a mild effect and translate the trend in /pair into a 
slightly more moderate mass dependence for the merger 
rate. 
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Fig. 5. — The major merger rate mass function in three redshift bins. The y-axis denotes ^, defined as a rate (per Gyr) per unit volume 
per logarithmic stellar mass interval. Solid circles indicate the observed merger rate for all galaxies determined with method I (see Fig[T]l 
and assuming the values of Tmg and Cmg from Patto n &l Atficld ( 2003) • Light blue circles show the result of excluding the approximate 
fraction of dry E/SO-E/SO mergers, as measured in ^ 15.31 The green triangles in the first redshift bin reflect the luminosity-dependent 
volumetric merger rate at z ~ 0.1 computed for 1:2 or greater m ass ratio mergers from the SDSS (Patton & Atfield 2008). Predictions 
for galaxy merger rates based on cosmological simulations from iStewart et al.l ||2008') are plotted as open diamonds. Red asterisk symbols 
show the rate of growth in the number of E/SO galaxies based on the sample from Bundy ct al. (2005). The rate of merging appears to 
be too low to fully account for the generation of new early-type galaxies. The dotted line in each panel indicates the average merger rate 
obtained if every galaxy experiences one merger over the range, 0.4 < z < 1.4. Rates lying in the shaded region above these lines have a 
strong impact on the population, while those below have a minimal effect. 



6.3. The Merger Rate and Merger Rate Mass Function 

Keeping in mind the substantial uncertainty in- 
volved with converting from /pair to a merger rate, 
we derive merger rates using the assumptions in 
both iKitzbichler fc White! (|2008l ) and iPatton fc Atfieldl 



()2008 ). In Table [3] we present the fractional galaxy 
merger rate, T^mg, defined as Tlmg, = Cmg/pair/Trng. T^mg 
can be thought of as the fraction of mergers per unit time 
within a samp le of galaxies. We take t he three most lu- 
minous bins in IPatton fc Atfieldl (|2008f l to correspond to 
our three mass bins and u se the fitting formula for 
(with Cmg set to 1.0) from IKitzbichler fc Whiti (|2008i ). 

With TZ-axg determined, we use our mass-limited sample 
to self-consistently compute the volumetric merger rate 
(the number of mergers per unit time per unit comoving 
volume) as a function of the stellar mass of the primary or 
host galaxy. We call this the merger rate mass function 
and denote it using the variable, ^ . The merger rate MF 
can be directly compared to differential evolution in the 
MFs of galaxies of various type, making it a valuable tool 



for evaluating the role of mergers in mass assembly and 
galaxy evolution. 

The merger rate mass function is determined by multi- 
plying the MF of the different host galaxy samples drawn 
from our survey (see 14. 4p by the relevant value of TZmg- 
In what follows, we use the values of TZmg determined 
from the merger timesc ales and efficiencies reported by 
IPatton fc Atfieldl (120081) and listed in Table El Using the 
values from IKitzbichler fc Whitd (|2008l ) leads to even 
higher merger rates (by ~10%) for the most massive 
galaxies. Our results are plotted as the dark solid cir- 
cles in Figure [5l The error bars reflect the uncertainty in 
/pair which dominates over the statistical uncertainty in 
the computed number densities. Both /pair and the cor- 
responding MFs have been determined using the same 
GOODS N-fS dataset described in § [21 It should be 
emphasized, however, that both suffer from cosmic vari- 
ance, although this problem is mitigated by combining 
two independent fields. We estimate that cosmic vari- 
ance affects the data plotted in Figure [5] at the 0.2 dex 
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level. The largest single source of uncertainty lies in the 
assumed values of Tmg and Cmg however. 

The volumetric merger rate of the SDSS sample (z 
0.1) measured in different l uminosity bin s (with pair 
separation 5-20 kpc) bv iPatton fc Atfieldl ()2008f ) is 
roughly plotted as a function of mass in the low-z bin 
of Figure [51 Given the lack of strong redshift evolution 
in our sample, the lower merger rate from SDSS may 
simply reflect the fact that IPatton k Atfieldl (2008) only 
consider merging pairs with an optical luminosity ratio 
greater than 1/2. Additional differences may arise in how 
Cmg and Tmg are applied to our different pair selection 
methods. Also shown in Figure [5] are predictions for the 
1:4 or greater M* ratio galaxy merger rate based on the 
cosmological N-body simulatio ns and halo modeling per- 
formed in St ewart et al.l ([lOOH)- These arc multiplied by 
the observed MFs to derive the volumetric rates plotted. 
The agreement is remarkable, especially given the un- 
certainties in Tmg. There is a hint that the merger rate 
derived from observations falls below the model predic- 
tions, especially at the highest masses. 

6.4. Hierarchical Mass Assembly Through Merging 

Figure [5] provides key insight about the nature of 
galaxy mass assembly. The dotted line in the figure 
provides a benchmark useful for gauging the impact of 
mergers on the galaxy population at different masses. 
It is determined by dividing the galaxy abunda nce at 
each redshift (from the MFs of iBundv et ani2006[ ) by 4.5 
Gyr, the amount of time spanning our redshift range, 
0.4 < z < 1.4. It is therefore the average event rate that 
would be obtained if every galaxy experienced exactly 
one merger over our redshift range. Rates that fall above 
this benchmark (in the shaded region of the plot) indi- 
cate a strong impact on the population. Rates that fall 
below correspond to processes with a minimal impact. 

Figure [5] shows that for galaxies with M^, < IO^^Mq, 
the observed merger rate lies below our benchmark, in- 
dicating that major mergers have not contributed sig- 
nificantly to the assembly history of such galaxies since 
z ^ 1. For more massive galaxies the situation changes. 
In the highest mass bin, the merger rate is roughly on 
par with the benchmark rate, demonstrating that major 
mergers play a larger role in the recent assembly of such 
galaxies. 

We can quantify this trend with the "merger remnant 
fraction," that is the fraction of systems that have under- 
gone major mergers in the different mass bins. This num- 
ber is determined by integrating TZ^g over 0.4 < z < 1.4. 
Using the results from method I, we find that at the 
highest masses, ~30% of galaxies experience a major 
merger. At the lowest masses probed, this number drops 
to 10-15%. These estimates are approximate because we 
have not considered transfe rs across mass bins as a re- 
sult of mergers and SF (see lDrorv fc Alvarejl2008f l and 
bec ause of the large uncertainties in Cmg and Tmg (also 
see lBell et al.l l2006bl. Still, they reinforce the increasing 
importance of mergers on the assembly history of higher 
mass galaxies at z ~ 1. 

Similarly, we can compare the volumetric merger 
rate to t he average SFR meas ured as a function of 
M, (e.g.. iDrorv fc Alvared l2008\ While measurements 
of SFR(M*) remain uncertain, this exercise indicates 
that the mass growth since z < 1.5 of galaxies with 



logM*/AfQ > 11 is almost completely dominated by 
merging. For systems with log M^/Mq ^ 10, how- 
ever, new growth from SF amounts to roughly 10 times 
the stellar mass accreted through major mergers, clearly 
indicating that SF is a far more important source of 
growth at lower masses. Considering the full mass range, 
logM*/MQ > 10, the mass accreted in major mergers 
over 0.4 < z < 1.4 amounts to ~15% of the total stel- 
lar mass density of systems with log M,,/Mq > 10 at 
z ~ 0.8. For log M^/Mq > 11, mass accretion through 
major mergers over the same redshift range accounts for 
-25% of at z - 0.8. 

Qualitatively, our observations of an increase in the 
merger fraction with mass echo the hierarchical assembly 
of dark matter halos as determined with ACDM N-body 
simulations. In a n ana lysis of the Millennium Simulation, 
iFakhouri fc Mai (|2008D . for example, find that the halo 
merger rate (for halo mass ratios greater than 1:3) rises 
by ~30% from Mdm = 10^^ to lO^^ This mass depen- 
dence is significantly weaker but in the same sense as the 
factor of —2 increase with mass we observe in the galaxy 
merger fraction among the 10 < logAf,/AfQ < 11.5 
galaxies thought to populate such halos. 

The enhanced mass dependence in the galaxy merger 
rate may be due to the way halos are occupied (details 
forthcoming in Hopkins et al., in prep), specifically the 
inferr ed peak in t he M^ /M^m ratio at MpM ~ lO^^Af© 
(e.g.. I Wang et all 120061 : iMandelbaum et al.ll2006t) which 
roughly corresponds to the low-mass end of the range 
probed by our sample. In this scenario, the major merger 
rate of halos and galaxies should be directly related near 
the M*/Mdm peak. But at higher masses Af* declines 
with respect to halo mass. This means that minor halo 
mergers can contribute to major galaxy mergers since less 
massive halos, being closer to the Mi,/M-£iyi peak, have 
larger Af* fractions. Because there are always more mi- 
nor than major halo mergers, this enhances the observed 
fraction of major galaxy mergers at higher masses. If our 
observations could probe below the peak, this interpre- 
tation would predict a suppressed galaxy major merger 
rate relative to the halo major merger rate since here 
some fraction of major halo mergers would host galaxies 
with a minor mass ratio. 

Indeed, this qualitati ve picture is borne out in detail 
bv lStewart et al.l ()2008D . As discussed above, the galaxy 
major merger rate predicted from their model based on 
N-body simulations is in remarkable agreement with our 
observations (see Figure [5]). 

6.5. Are Spheroidals Formed in Major Mergers? 

Another important insight gained from Figure [5] is the 
role of mergers in the formation of elliptical and red 
sequence galaxies. The red asterisks in this plot show 
the formation rate (the rate at which the abundance in- 
creases) of visually classified spheroidal sys t ems, as de- 
termined from the sample of iBundv et al.l ()2005f ). Er- 
ror bars are Poissonian and do not include cosmic vari- 
ance which we estim ate enters at the —0.2 dex level (see 
iStringer et aTll200l . 

Figure |S| shows that the formation rate of new 
spheroidal galaxies is greater than the merger rate at 
nearly all masses and redshifts probed. The disagree- 
ment increases when only wet or mixed mergers are 
considered — only these can create new spheroidals. The 
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hght blue points show the merger rate mass function af- 
ter the approximate fraction of dry spheroidal-spheroidal 
mergers (as measured in § 15. 3p is subtracted. This re- 
sult was anticipated by an analysis of the expected num- 
ber of major halo mergers as measured in the Millen- 
nium Simulation which is also unable to account for the 
rate at which new spheroidals appear (jBundv et al.ll2007l : 
iGenel et al.ll2008l ) . A similar conclusion applies to the in- 
ability of major mergers to account for the formation of 
red sequence galaxies which increase at a similar rate as 
the spheroidal population fe.g.. ,Bell et al.ii2007i' ). 

The discrepancy between major (wet) mergers and the 
formation of spheroidals shown in Figure [5] ranges from 
factors of 3-12. Simply interpreted, this requires 
mechanisms in addition to merging to make spheroids. 
How strong is this requirement? One concern is cos- 
mic variance which could affect both the merger and 
spheroidal formation rates — we estimate the discrepancy 
could be overestimated by no more than a factor of ~2 as 
a result (Strin ger et al.l[2008f l. At the same time, the gap 
could be closed to some extent if the merger timescale 
were less than wO.5 Gyr (e.g., iLotz et al] 1200 8a'). In- 
deed, a comparison of various theoretical treatments con- 
cludes that current predictions for T^g are again uncer- 
tain at the factor of ^2 level (Hopkins et al., in prep). 
Furthermore, our spheroidal classification, while domi- 
nated by ellipticals, includes SO's and some Sa galax- 
ies^^. The large bulges in such galaxies could be built 
in mergers with mass ratios less than 1:4, increasing the 
number of relevant mergers by 30%. Finally, multiple 
minor mergers can also drive morphological evolution. 
We would expect a factor of ^2 times more 1:10 minor 
mergers than the nu mber of observed 1:4 mergers (e.g., 
I Stewart et al.l l2008f ) . If ^--^3 such minor mergers have 
the equivalent morphological impact as one 1:4 merger 
(jBournaud et al. 2007; Hopkins et al. 2008), minor merg- 
ers can account for an additional 2/3 « 70% of the mor- 
phological transformations needed to form spheroidals. 

Assuming these potential systematic effects conspire 
in the same way and including minor mergers, it would 
be possible to reduce the discrepancy by a factor of ~6. 
This would lend stronger support to the merger hypoth- 
esis and motivate a more precise treatment, but would 
still leave some tension at z ~ 0.8. For this reason and 
given the current observations and most recent theoret- 
ical assumptions, we conclude that other mechanisms 
in addition to merging may be needed to explain the 
transformation of disk-like galaxies into spheroidals (see 
iParrv et al.ll200^ and the quenching of star formation. 

7. SUMMARY 

Using a deep i^g-selected catalog comprised of data 
from the GOODS fields, we have presented an analysis of 
the close pair fraction and implied merger rate of galax- 
ies drawn from a mass-limited sample over the redshift 
range 0.4 < 2: < 1.4. Using two methods of estimation, 
we find a relatively low pair fraction of '^4% with a red- 
shift dependence of (1 -f- z)^-^^^-^. Our analysis strongly 
supports the basic conclusion of a low merger rate since 
z 1, as deduced by several other recent surveys. 

The fact that it is so difficult to discriminate between such 
galaxies at 2 > 0.4 i s the reason we group t hem into a single 
"spheroidal" category IIBrinchmann et aLlliggsf) . 



Our key observational finding is that the pair fraction 
of galaxies with Af* > lO-^^M© is higher (/pair ~ 5-9%) 
than the corresponding pair fraction of lower mass sys- 
tems (/pair ~ 2-4% for ~ IO^^Mq). In addition, red 
systems and galaxies with spheroidal morphologies ap- 
pear more likely to have companions than their blue or 
disk-like counterparts. This is in line with extrapolations 
of the correlation function to small scales. We find that 
the fraction of host galaxies classified as irregular (often 
more than half of hosts with Af* < 3 x 10^°) is larger 
than the fraction of irregulars in the parent population 
(typically 10-20%). This supports the notion that or- 
biting companions can distort the primary galaxy even 
before a merger has occurred. 

Dry mergers are only significant at the highest masses 
in our sample, > IO^^Mq, but this appears to result 
from the fact that most galaxies in this mass range have 
already become red, early- type systems by z ~ 1. At 
lower masses (Af* ^ lO^^M©) the fraction of red galaxies 
that host companions is consistent with zero, and we thus 
conclude that dry mergers alone are incapable of fully 
explaining the build-up of the red sequence, especially at 
M^, « 2-3 X IO^^Mq where the formation of red galaxies 
is particularly rapid. In support of this conclusion, we 
also find that the number of dry pairs accounts for ~50% 
of all pairs in our mass range and evolves in way that is 
consistent with pairs drawn randomly from the parent 
population. The emerging picture is one in which dry 
mergers are not the main driver for red sequence growth 
but instead become increasingly common as a result of 
it. 

By adopting estimates for the close pair merger effi- 
ciency and timescale, we determine the mass-dependent 
volumetric merger rate. From a comparison to the for- 
mation rate of galaxies with spheroidal morphology at 
the same redshifts, we conclude that the major merger 
rate is too low to fully explain the formation of spheroidal 
(and red sequence) galaxies since z ~ 1. 

In terms of mass assembly, major mergers have a 
strong impact on galaxies with Af* > IO^^Mq, ~30% 
of which have experienced a major merger during the 
time spanning our redshift range. For less massive galax- 
ies (AT* < 3 X IO^^Mq), such mergers play a less sig- 
nificant role, affecting less than 10% of the population. 
For systems at such intermediate and lower masses, star 
formation is likely to be of much greater importance in 
driving growth. The higher merger fraction observed for 
the more massive galaxies in our sample provides direct 
support for hierarchical assembly and agrees with recent 
galaxy models based on N-body simulations, even if ad- 
ditional processes are required to explain the top-down 
nature of star formation downsizing and the formation 
of spheroidal galaxies. 
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APPENDIX 
MOIRCS REDUCTIONS 

We provide additional information on the MOIRCS data reduction, which was carried out using a modified version 
of the MCSRED package written by Ichi Tanaka. Simply using the package as written often led to the appearance 
of fringes across the entire reduced frame, especially for detector 2. More than half of our pointings were strongly 
affected, and the typical pcak-to-trough fringe height was a few percent of the background level. 

We traced the cause to two problems. First, strong fringing in one or two frames in an image set would be imprinted 
into other frames via the skyflat procedure and second, skyflats in MCSRED are made by averaging all images in a 
set. The flat field can vary over timescales of minutes, however, and changes (as well as the fringing) may be related to 
dithering the telescope. Thus a "running" sky flat that includes only some images taken before and after the science 
frame is desirable. We wrote a separate procedure which constructs running skyflats to flatfield the data and allows 
frames with bad fringing to be removed. The image registration portion of MCSRED also caused occasional problems, 
with some images being excluded from the final mosaic. By default, MCSRED calls the IRAF routine xyxymatch 
with the automated "triangles" registration option, which can sometimes fail for certain frames. When this happened 
we identified the positions of three stars common to all frames in an image set and supplied these to xyxymatch 
using the "tie points" method. 

Careful inspection of the MOIRCS images revealed a curious pattern of noise spikes (although "noise worms" is a 
better description since they tend to be elongated) that often appear within ~10" of bright sources (often stars) with 
Kg < 16. These objects are often detected as sources with Kg « 23 and, when they appear near galaxies, can be 
mistaken for fainter companions. We identified 133 such noise worms (about 1% of MOIRCS detections) and removed 
them from the photometric catalogs. 
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TABLE 1 

ReSFLTS for /p„;i- 



logM^/M0 > 10 10 < log M«/M0 < 10.5 10.5 < logM,/MQ < 11 logM^/Mg > 11 



Sample 


z 


Np 


^corr 


Nqal 


/pair (%) 


Np 


^corr 




/pair (%) 


Np 


^corr 




/pair (%) 


Np 


^corr 


Nqal 


/pair (%) 


Method I 


AH 


0.4-0.7 


60 


43.9 


514 


3 it 2 


37 


29.4 


250 


3 ± 2 


18 


12.7 


200 


3 ± 2 


5 


1.8 


64 


5 ± 2 


All 


0.7-0.9 


56 


36.5 


380 


5 ± 3 


33 


23.6 


183 


5 ± 3 


16 


10.4 


136 


4 ± 3 


7 


2.6 


61 


7 ± 3 


All 


0.9-1.4 


126 


79.8 


793 


6 ± 2 


53 


41.7 


326 


3 ± 2 


54 


30.9 


333 


7 ± 2 


19 


7.1 


133 


9 ± 2 


Blue 


0.4-0.7 


33 


24.5 


268 


3 ± 3 


27 


20.0 


183 


4 ± 3 


6 


4.2 


75 


2 ± 3 





0.3 


10 


-3 ± 3 


Blue 


0.7-0.9 


33 


21.9 


212 


5 ± 3 


28 


17.4 


142 


7 ± 3 


4 


4.0 


58 


± 3 


1 


0.4 


12 


5 ± 3 


Blue 


0.9-1.4 


87 


59.2 


574 


5 ± 2 


47 


38.4 


304 


3 ± 2 


33 


17.8 


206 


7 ± 2 


7 


3.1 


63 


6 ± 2 


Red 


0.4-0.7 


27 


19.4 


246 


3 ± 3 


10 


9.4 


67 


1 ± 3 


12 


8.5 


125 


3 ± 3 


5 


1.5 


54 


6 ± 3 


Red 


0.7-0.9 


23 


14.6 


168 


5 ± 4 


5 


6.2 


41 


-3 ± 4 


12 


6.3 


78 


7 ± 4 


6 


2.2 


49 


8 ± 4 


Red 


0.9-1.4 


39 


20.5 


219 


8 ± 4 


6 


3.3 


22 


12 ± 4 


21 


13.1 


127 


6 ± 4 


12 


4.1 


70 


11 ± 4 


Method II 


All 


0.4-0.7 


21 


8.2 


514 


3 ± 1 


8 


5.5 


250 


2 ± 1 


10 


2.3 


200 


5 ± 1 


3 


0.3 


64 


5 ± 1 


All 


0.7-0.9 


23 


10.1 


380 


4 ± 2 


10 


6.6 


183 


3 ± 2 


10 


2.9 


136 


6 ± 2 


3 


0.6 


61 


5 ± 2 


All 


0.9-1.4 


45 


22.5 


793 


4 ± 1 


15 


10.9 


326 


2 ± 1 


20 


9.6 


333 


4 ± 1 


10 


2.0 


133 


7 ± 1 



Note. — Wc define Np as the raw number of pairs, while Ncorr is the estimated number of contaminants. The number of host galaxies is given by Ngai- Method II /| 
values have been corrected upward by 1% to account for catastrophic photo-2 errors. 



